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Abstract 

Background: Serum response factor (SRF) is a widely expressed transcription factor involved in nnultiple regulatory 
programs. It is believed that SRF can toggle between disparate programs of gene expression through association 
with different cofactors. However, the direct evidence as to how these factors function on a genome-wide level is 
still lacking. 

Results: In the present study, I explored the functions of SRF and its representative cofactors, megakaryoblastic 
leukemia 1/2 (l\/lKLl/2) and ETS-domain protein 4 (ELK4), during fungal infection challenge in macrophages. The 
knockdown study, combined with gene expression array analysis, revealed that )\/lKLl/2 regulated SRF-dependent 
genes were related to actin cytoskeleton organization, while ELK4 regulated SRF-dependent genes were related to 
external stimulus responses. Subsequent chromatin immunoprecipitation coupled with massively parallel sequencing 
(ChlP-seq) suggested that many of these regulations were mediated directly in cis. 

Conclusions: I conclude that SRF utilizes l\/lKLl/2 to fulfill steady state cellular functions, including cytoskeletal 
organization, and utilizes ELK4 to facilitate acute responses to external infection. Together, these findings indicate that SRF, 
along with its two cofactors, are important players in both cellular homeostasis and stress responses in macrophages. 

Keywords: SRF, MKLl, )\/lKL2, ELK4, Gene expression array, ChlP-seq 



Background 

Macrophages are an important subset of innate immune 
cells implicated in host defense against bacterial and viral 
infections. Three of their major functions include phago- 
cytosis and degradation of foreign pathogens, clearance of 
apoptotic or necrotic host cells, and rapid generation of 
effector cytokines that orchestrate immune responses. 
These different functions of macrophages are likely regu- 
lated by distinct transcription programs. 

Serum response factor (SRF) is a member of the MADS 
(Mcml, Agamous, Deficiens, 5'RF) transcription factor 
family [1] and is abundantly expressed in a variety of cell 
types. Deletion of Srf is embryonic lethal [2]. Conditional 
deletion of Srf revealed its essential functions in many 
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tissues, including skeletal, cardiac and smooth muscle cells 
[3,4], neurons [5], and thymocytes [6,7]. At the cellular 
level, these phenotypes can be attributed to impaired 
expression of different sets of SRF target genes. SRF- 
regulated immediate early genes (lEGs), such as c-fos, 
early growth response 1 {Egrl) and junB, are important 
for cell growth and proliferation [8-10], while SRF- 
regulated troponin, vinculin (Vc/), and SM22a partici- 
pate in muscle-specific contractile functions [9,11]. 

With high affinity and specificity, homodimeric SRF 
binds to a consensus c/5-element known as a CArG box, 
with the sequence of CC{A/T)^GG [12,13]. SRF regulates 
diverse programs of gene expression through its associ- 
ation with different accessory factors. The most well- 
studied cofactors of SRF include members of the ternary 
complex factor (TCF) and myocardin families [9] of tran- 
scription factors. The TCF family of proteins includes 
ELKl, ELK3 (also known as SAP2, Net) and ELK4 (also 
known as SAPl). TCFs respond to mitogen-activated pro- 
tein kinase (MAPK) signaling initiated by serum or growth 
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factors [14] and is involved in the regulation of lEGs. 
TCFs interact with SRF and simultaneously bind to an 
ETS binding site adjacent to the CArG box. Together, the 
CArG box and the ETS site are recognized as the serum 
response element (SRE) [15]. Myocardin family members 
include myocardin, megakaryoblastic leukemia 1 (MKLl, 
also known as MAL, MRTF-A and BSAC) and MKL2 
(also known as MRTF-B). Myocardin is strictly expressed 
in cardiac and smooth muscle cells, while MKLl and 2 
exhibit broader expression patterns [16]. MKLs connect 
actin dynamics with SRF activation [17] and can regulate 
muscle-specific contractile genes in response to actin 
polymerization [18]. Unlike TCFs, MKLl and MKL2 do 
not bind DNA sequences directly [19]; instead, they associ- 
ate with SRF to regulate the transcription of target genes. 
SRF interaction with MKLs and TCFs is mutually exclusive 
[20,21], but recently, MKLs were reported to regulate the 
expressions of some lEGs in response to serum induction 
[22,23]. However, little is known about how these two fam- 
ilies coordinate with SRF on a genome-wide level. 

Macrophages express a variety of toll-like receptors 
(TLRs) on their cell surface that enable recognition of 
molecules, known as pathogen-associated molecular pat- 
terns (PAMPs) [24], that are broadly expressed by patho- 
gens, but not by the host. Zymosan, a yeast cell wall 
derivative, induces inflammatory responses in macro- 
phages through binding to the dectin-1 receptor or TLR2/ 
6 heterodimers [25]. The function of SRF in macrophages 
had not been investigated until our previous studies 
revealed that SRF regulates both general and cell type- 
restricted cytoskeletal gene expression programs in these 
cells [26]. Interestingly, we also found that the SRF/MADS 
(CArG box) motif is highly enriched in the promoters of 
zymosan-induced genes and that SRF indirectly modulates 
type I interferon signaling in macrophages [27]. Overall, 
these studies suggest a possible role of SRF in the immune 
responses of macrophages. 

With recent advances in paraUel high throughput se- 
quencing, I sought to map out the genome-wide occu- 
pancy of SRF and its cofactors, MKLl/2 and ELK4. In 
this study, RNA interference was used to knock down 
SRF, MKLl/2, and ELK4, foUowed by genome-wide ex- 
pression profiling analysis and ChlP-seq experiments to 
identify the target genes and direct binding sites of these 
cofactors. Together, I demonstrate that the two SRF co- 
factor families, TCFs and MKLs, regulate distinct subsets 
of SRF target gene programs during fungal infection in 
macrophages. 

Results 

SRF is required for macrophage responses to 
zymosan treatment 

In order to identify the transcription factors that are in- 
volved in gene regulation in response to fungal infection 



in macrophages, I stimulated primary macrophages with 
zymosan for 1 h and harvested RNA for microarray ana- 
lysis. Motif analysis on the promoters of genes shown by 
the array to be induced by zymosan identified the SRF/ 
MADS (CArG box) motif as one of the most highly 
enriched motifs [27] (Additional file 1: Figure SI A). This 
is in accordance with the fact that SRF functions as an 
immediate responder to external stimuli such as growth 
factors or serum. Interestingly, Srf transcription was 
increased 2-fold after 1 h of zymosan treatment, but 
returned to baseline after 6 h of treatment (Additional 
file 1: Figure SIB), suggesting that SRF activation is an 
immediate early phase response. To evaluate if this 
change in mRNA expression had a functional effect, I 
transfected RAW264.7 cells with a reporter construct 
containing 3 tandem SRF-binding c/5-elements (SxCArG). 
After 1 h of treatment with zymosan, I found an approxi- 
mately three-fold increase in reporter activity (Additional 
file 1: Figure SIC). Next, I evaluated the impact of knock- 
ing down Srf in zymosan treated elicited mouse peritoneal 
macrophages by using microarray analysis [27]. Consider- 
ing that the transcription of Srfwzs, induced quickly and 
that the SRF motif was enriched in the early phase of 
zymosan treatment, I decided to focus the rest of my 
study on cells treated with ligand for short periods of time. 
To identif)^ SRF-dependent genes, I filtered the microarray 
results to isolate the set of transcripts whose expressions 
were reduced by more than 1.5 fold with SRF-specific 
siRNA after 1 h of zymosan treatment. The resulting sub- 
set contained 264 non-redundant downregulated genes 
and 203 upregulated genes (full list in Additional file 2: 
Table SI), which indicates that SRF can be induced by 
zymosan stimulation and that it is an important contribu- 
tor to macrophage responses to external stimuli. 

To investigate which biological functions the SRF target 
genes are involved in, I performed gene ontology (GO) 
analysis of the SRF regulated genes using the DAVID 
Functional Annotation Tool. The top 5 functional annota- 
tions are shown in Figure lA (fuU results in Additional 
file 3: Table S2, Part I). The most significantly enriched 
item was "actin cytoskeleton organization", which is a 
function of SRF that has been elucidated in other cell 
types and in our previous study in quiescent macro- 
phages [26] . It is noteworthy that some genes within this 
subgroup , such as Corola [28], Lspl [29] andLcpi [30], 
are preferentially expressed in cells that are of myeloid 
or hematopoietic origin and regulate macrophage specific 
functions such as phagocytosis or migration. Moreover, 
other GO terms that are closely related to macrophage 
specific functions were also significantly enriched such as 
"regulation of cell adhesion" (including Cd24a, Cited2, 
Csfl, Itga6, MMP14, Srf, Thbsl and Tgm2) and "response 
to external stimulus" (including Cd24a, BatS, Hpsl, MEFV, 
Argl F13AI CIQB, CFP, Corola, ENTPDl, GPR77, Igfbp4, 
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Figure 1 SRF is required for macrophiage function in response to zymosan stimulation. (A) Gene ontology analysis for biological process 
annotations of downregulated and upregulated SRF dependent genes after zymosan treatment for 1 h. (B) mRNA expression profiles of the 
indicated representative SRF dependent genes were analyzed by Q-PCR. Thioglycollate-elicited macrophages were transfected with nonspecific 
(NS) or specific siRNAs for SRF and treated with or without zymosan for 1 h. """p < 0.01 vs. NS; ""p < 0.05 vs. NS. Target gene expression is 
represented as the average fold change compared to NS siRNA transfected samples of at least three independent experiments. Error bars 
represent standard deviations. 



Igta6, Lcpl Lspl MmplS, Mmpl4, Mest, USF2, Thbsl 
and IllO), 

I validated the expression changes of some representa- 
tive target genes by quantitative PGR (Q-PCR). As ex- 
pected, several groups of target genes were confirmed as 
SRF dependent, including: (1) general actin cytoskeleton 
related genes such as cofilinl {Cfll) and Vcl, and the mi- 
gration related gene, lymphocyte-specific protein 1 {Lspl); 

(2) adhesion related genes like thrombospondin 1 {Thbsl); 

(3) genes related to stimulus responses such as IllO and 
EgrS; and (4) genes related to other functions, including 
/m5-related tyrosine kinase 1 {Fltl) and fos-like antigen 2 
{Fosl2) (Figure IB). All of these results indicate that SRF is 
required for intact expression of genes closely related to 
macrophage functions. 

GO analysis of the genes upregulated by SRF knockdown 
under zymosan treated conditions resulted in the identifica- 
tion of more generic terms such as "cellular process", 
"nucleosome assembly" and "regulation of metabolic process" 
(Figure lA, full results in Additional file 3: Table S2, 
Part II). Besides regulating the transcription of target 
genes, SRF is known to regulate some miRNAs [31-33] 
associated with altered expression of mRNAs and pro- 
teins. The genes upregulated upon SRF knockdown could 
potentially be downstream targets of SRF through the 
regulation of miRNAs. Using Ingenuity Pathway Analysis, 



I compiled a list of 10 overrepresented miRNAs for 
the 203 upregulated genes, which included miR-16-5p, 
miR-489, miR-205-5p, miR-101, miR-291a-3p, miR-204- 
5p, miR-33, miR-615-3p, miR-424-3p and miR-10a-5p. 
Although there is no current literature linking any of 
these miRNAs to SRF, it would be interesting to test, in 
future studies, whether these miRNAs are actually novel 
targets of SRF. 

MKL1/2 and ELK4 regulate different subsets of genes in 
zymosan treated macrophages 

Previous studies of SRF in muscle cells suggested that SRF 
participates in different regulatory programs by associat- 
ing with different transcription cofactors. This mechanism 
enables SRF to toggle between regulating cell proliferation 
and growth and orchestrating actin cytoskeletal and con- 
tractile homeostasis. I next sought to investigate accessary 
factors utilized for SRF function in response to zymosan 
stimulation in macrophages. I focused on two well-known 
SRF associated cofactor families: the myocardin and TCF 
families of transcription cofactors. 

MKLl and MKL2, share many redundant functions [23] 
and are both highly expressed in primary macrophages. 
Unlike Srfl mRNA levels of Mkll and Mkl2 stay constant 
after zymosan treatment (Additional file 4: Figure S2A). 
Because of their many redundancies, I used an siRNA 
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pool targeting both Midi and Mkl2 to knock down their 
expressions simultaneously to better assess their role in 
macrophages (Additional file 4: Figure S2A, MKLl/2 will 
denote both MKLl and MKL2 in the following text). Elk4 
is the dominantly expressed member of the TCF family 
in primary macrophages (Additional file 4: Figure S2B). 
Similar to Mkll/2, Elk4 expression does not change dur- 
ing zymosan treatment (Additional file 4: Figure S2B). 
Transfection of macrophages with H/c^-specific siRNAs 
significantly reduced Elk4 expression compared with 
control siRNA transfected samples (Additional file 4: 
Figure S2C). I also confirmed that the knockdown effect 
is specific for Elk4 and that mRNA expression for Elkl 
and Elks were unaffected (Additional file 4: Figure S2C). 
Following confirmation of the siRNA knockdown effi- 
ciencies, I performed gene expression array analysis on 
my samples to compare the gene expression changes in 
Mkll/2 and Elk4 knockdown cells versus control siRNA 
transfected cells after zymosan treatment for 1 h. 

Using 1.5 fold reduction as a cutoff, I identified a total 
of 309 genes that were MKLl/2 dependent (full list in 
Additional file 5: Table S3) and 265 genes that were ELK4 
dependent (full list in Additional file 6: Table S4). These 



A 



results indicate that the absence of these cofactors would 
have an impact on zymosan induced gene expression in 
macrophages. 

To check if MKLl/2 contributes to SRF related functions 
in macrophages, I compared the 309 MKLl/2 target genes 
with the 264 SRF target genes. This analysis identified 60 
MKLl/2 target genes that were also dependent on SRF for 
their expression (Figure 2A). GO analysis of these 60 genes 
identified the term of "actin cytoskeleton organization" as 
the top hit. When I performed GO analysis on the MKLl/2 
target genes that were independent of SRF expression, I 
found that "positive regulation of biological process" ranked 
as the first significant term (Figure 2B). Tal<en together, 
these results suggest that MKLl/2 could coordinate with 
SRF to regulate cytoskeleton related target genes, but may 
also have other generic functions independent of SRF. 

By comparing the 265 ELK4 target genes with the 264 
SRF target genes, I identified 24 genes that were dependent 
on both ELK4 and SRF for their expression (Figure 2A). 
GO analysis of these 24 genes identified two significantly 
enriched terms: "response to external stimulus" and "regu- 
lation of transcription from RNA polymerase II promoter". 
This result indicates that ELK4 is involved in a different 
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Figure 2 MKLl/2 and ELK4 regulate the transcription of different gene sets in zymosan induced macrophages. (A) Venn diagram 
indicating tine relationsliips between SRF dependent genes witli respect to MKLl/2 and ELK4 dependence. (B) Gene ontology analysis for 
biological process annotations of the SRF dependent or independent MKLl/2 target genes after zymosan treatment for 1 h. (C) Gene ontology 
analysis for biological process annotations of the SRF dependent or independent ELK4 target genes after zymosan treatment for 1 h. 
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aspect of SRF function than MKLl/2. The fact that ELK4 
dependent genes are enriched for the GO term related to 
stress responses in gene ontology analysis is consistent with 
our finding that besides the SRF motif, the ELK4 motif is 
also highly enriched in the promoters of zymosan-induced 
genes in primary macrophages (Additional file 1: Figure 
SIA). Since the majority of ELK4 dependent target genes 
are not dependent on SRF for expression, ELK4 may have 
some SRF independent functions in macrophages. Simi- 
lar GO analysis on the SRF independent ELK4 target 
genes showed the term "cellular process" shifted to the 
top, followed by "response to stress" and "cell division", 
etc. (Figure 2C). 

Taken together, these findings support the view that 
MKLl/2 and ELK4 are involved in distinct functions of 
SRF in zymosan treated macrophages. Specifically, MKLl/2 
is involved in actin cytoskeleton organization, ELK4 con- 
tributes to establishing responses to external stresses, and 
both sides constitute important functions of SRF in macro- 
phages. Moreover, I found that both MKLl/2 and ELK4 



have roles that are independent of SRF and that they may 
perform these roles independently or in coordination with 
other transcription factors. Additional work will be ne- 
cessary to explore how SRF and these cofactors work on 
the genome directly and also to identify other potential 
cofactors. 

MKL1/2 and ELK4 dependent target genes show different 
responses to zymosan induction 

Similar to the SRF target genes, I also validated repre- 
sentative target genes of MKLl/2 and ELK4 by Q-PCR. I 
verified the MKLl/2 dependence of a subset of genes 
with functional annotations linked to actin cytoskeleton 
organization, including Cfll, Dstn, Vcl, and Thbsl, All four 
of these genes are highly expressed in macrophages in the 
basal state. While the former two genes were not regu- 
lated by zymosan, the latter two were slightly upregulated 
by zymosan treatment. Expression of all 4 genes was re- 
duced 30-60% by MKLl/2 knockdown (Figure 3A). More- 
over, ELK4 expression was not required for the expression 
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Figure 3 Expression profiles of l\/lKL1/2 and ELK4 dependent target genes in response to zymosan treatment. (A) and (B) Expressions of 
MKLl/2 dependent genes {Cfll Dstn, Vcl and Thbsl) and ELK4 dependent genes {IllO, c-fos, Egrl and EgrS) were analyzed by Q-PCR. Thioglycollate-elicited 
macrophages were transfected with NS or specific siRNAs for MKLl/2 or ELK4, as indicated, and treated with or without zymosan for 1 h. (C) Srf mRNA 
expression in response to MKLl/2 or ELK4 l<nocl<down under basal conditions or after treatment with zymosan for 1 h. Target gene expression is 
represented as the average fold change compared to NS siRNA transfected samples of three independent experiments. **p < 0.01 vs. NS; *p < 0.05 vs. NS. 
Error bars represent standard deviations. 
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of genes related to actin cytoskeleton organization such as 
Vcl and Cnn2 (Additional file 4: Figure S2D). 

I also verified the ELK4 dependence of a subset of stimu- 
lus response genes, including Egrl, c-fos, IllO, and EgrS 
(Figure 3B). The expression of these four genes can be in- 
duced as much as 100 fold by zymosan, in contrast to actin 
related genes. Although not revealed by microarray evalu- 
ation, I found by Q-PCR that knocking down MKLl/2 also 
impaired the expression of some stimulus response genes 
such as EgrS and II 10 (Additional file 4: Figure S2E). 

Generation of a BLRP-MKL1-BirA stable cell line for 
ChlP-seq analysis in RAW264.7 macrophages 

While knockdown experiments combined with mRNA 
expression analysis helped to identify the genes whose 
expressions are dependent on SRF, MKLl/2 and ELK4, it 
could not tell me whether these dependencies are through 
direct transcriptional regulation or through a secondary 
effect. It has been reported (and our expression data con- 
firms) that SRF regulates the expression of itself [34,35] 
and that this effect is dependent on both MKL and TCF 
cofactors [14,18,36]. Thus, it is possible that knocking 
down MKLl/2 and ELK4 could cause reduced expressions 
of some target genes through reducing the expression of 
Srf, As shown in Figure 3C, knocking down MKLl/2 
reduced the expression of Srf more than 50%, under 
both basal and zymosan treated conditions (p < 0.001 and 
p = 0.017, respectively). Knocking down ELK4 showed a 
similar effect (p = 0.036 and 0.013, respectively). This re- 
sult is consistent with previous reports that there is an 
MKL and TCF-dependent positive feedback loop for 
SRF-dependent autoregulation [14,18,36]. 

To identify the direct targets of SRF and its cofactors on 
a genome-wide level, I performed ChIP followed by high 
throughput sequencing of the enriched DNA fragments 
(ChlP-seq). MKLl and MKL2 do not bind DNA sequences 
directly [19]; instead, they associate with SRF to regulate 
the transcription of target genes. As the binding affinity of 
MKLl/2 to the DNA fragment is relatively weak, conven- 
tional Chip assay with MKLl antibody failed to provide 
sufficient enrichment on known target genes in macro- 
phages (data not shown). To this end, I generated a BLRP- 
tagged MKLl expression vector and transfected it into 
RAW264.7 cells to generate a stable cell line exoge- 
nously expressing BLRP-MKLl (BLRP-MKLl-BirA 
stable cell line). A TEV-based ChIP strategy was then 
performed for MKLl using this stable cell line, followed 
by sequencing analysis. I note here that SRF and ELK4 
directly bind DNA, so ChlP-seq using commercial anti- 
bodies against endogenous proteins were possible in the 
MKLl overexpressed BirA stable cell line. 

Stable clones with highly expressed BLRP-MKLl were 
identified and selected by RT-qPCR with primers targeting 
the BLRP sequence (Additional file 7: Figure S3 A). Blrp 



expression was approximately 4000 times higher in the 
stable cell line compared to parental BirA cells and the 
Blrp mRNA can be further induced as much as 10 fold 
following 1 h of zymosan treatment. MKLl expression 
levels in this stable clone were twice as high as in the BirA 
parental cell line under basal conditions and were in- 
duced 7.5 fold after 1 h of zymosan treatment (Additional 
file 7: Figure S3A). The high expression of MKLl protein 
was also confirmed by western blot (Additional file 7: 
Figure S3B). 

Global localization of SRF, MKL1 and ELK4 in macrophages 

All the ChlP-seq assays were performed in zymosan treated 
BLRP-MKLl-BirA stable cells. Using a false discovery rate 
of 0.1%, I identified a total of 5025 SRF-specific peaks, 1140 
MKLl -specific peaks and 783 ELK4-specific peaks (Full 
lists are provided in Additional file 8: Table S5, Additional 
file 9: Table S6 and Additional file 10: Table S7, respec- 
tively). Each of these peaks were then annotated to the clo- 
sest RefSeq gene observed. Among the 5025 SRF-bound 
peaks, I found 781 and 350 peaks that were also associated 
with MKLl and ELK4 binding, respectively. Furthermore, 
104 SRF peaks were associated with both MKLl and ELK4 
(Figure 4A). There were also a number of SRF associated 
peaks that did not overlap with MKLl and ELK4 binding 
sites, indicating that there may be other cofactors utilized 
by SRF. For MKLl, the majority of the peaks (68.5%, 781 
out of 1140) were associated with SRF. However, the per- 
centage for ELK4 association with SRF peaks was relatively 
lower (44.7%, 350 out of 783), indicating that ELK4 may 
have other functions independent of SRF and could have 
other co-regulators [37]. 

Motif analysis of the 5025 nearest genes associated 
with SRF peaks by oPOSSUM recovered many tran- 
scription factor binding sites (TFBSs). As expected, SRF 
and ELK4 were among the top hits, together with GABPA, 
NFkB, Staf, TEADl, NFYA, RELA, CREBl, STATl and 
RREBl (Figure 4B). Some of the overrepresented TFs have 
been reported as cofactors of SRF already, including NFkB, 
TEADl, NFYA and RELA [38-42]. My study also revealed 
additional potential cofactors of SRF, such as GABPA, 
Staf, CREBl, STATl and RREBl. The most highly en- 
riched motif for MKLl bound genes was SRF, followed by 
MEF2A, NFkB, STATl, IRFl, MIZF, FOXF2, GABPA, 
NFYA, RELA and CREBl (Figure 4B). Most of the over- 
represented TFBSs of MKLl peaks coincided with those 
of SRF (including NFkB, STATl, GABPA, NFYA, RELA 
and CREBl), suggesting that MKLl might be involved in 
the same complex with SRF and other cofactors. However, 
there is also the possibility that MKLl might interact with 
the novel cofactors, independently of SRF, to regulate spe- 
cific target genes. 

For ELK4 bound genes, ELK4 and SRF motifs were 
significantly enriched, together with NFYA and GABPA 
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Figure 4 Genomic analysis of SRF, MKLl and ELK4 ClilP-seq peaks in macropiiages. (A) Venn diagram of the overlap between SRF, MKLl, 
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(Figure 4B). NFYA has been predicted to interact with 
both ELK4 and ELKl [43]. GABPA, another ETS-domain 
transcription factor, has been reported to show a substan- 
tial overlap in promoter occupancy with ELKl and ELK4 
[44]. The association of ELK4 with both NFYA and 
GABPA motifs suggests that, besides SRF, these factors 
could be new partners of ELK4 in macrophages and 
should be experimentally validated in future work. 



Classification of SRF peaks with respect to the loca- 
tions of genes revealed that a large fraction of SRF peaks 
were located in intergenic regions (37.1%) and introns 
(32.4%). Approximately 27.4% of the peaks were located 
in the vicinity of promoters (+/- 1000 bp of a transcrip- 
tion start site) and only 3.0% in exons (Figure 4C). This 
distribution pattern is consistent with our previous SRF 
ChlP-seq studies in quiescent primary macrophages [26]. 
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A large fraction of the peaks seen in primary macrophages 
(923 peaks, 73.1% of the total SRF bound peaks) were also 
identified in my SRF ChlP-seq in BLRP-MKLl-BirA cells. 
Genomic annotation of MKLl binding sites indicated a 
distribution pattern similar to SRF peaks. 39.7% of the 
MKLl peaks were located in intergenic regions, 30.5% in 
introns, 26.8% in promoters and a very small fraction in 
exons (3.0%) (Figure 4C). However, the genomic distribu- 
tion of ELK4 binding sites differed significantly from that 
of SRF and MKLl. Approximately 42.4% of the ELK4 
peaks were located in promoters, 32.1% in intergenic re- 
gions, 21.2% in introns and a very small fraction in exons 
(4.4%) (Figure 4C). This data indicates that ELK4, unlike 
SRF and MKLl, may function primarily through promoter- 
proximal regions. 

The fact that there is a difference in the number of 
SRF ChlP-seq peaks in intergenic regions versus those of 
ELK4, indicates that SRF participates in some unique 
transcription regulation mechanisms independent of its 
well-known cofactor, ELK4. There are two possible expla- 
nations for this phenomenon. First, SRF could function at 
distal enhancers for its regulated genes. This possibility 
was verified by our previous study showing that SRF binds 
together with PU.l at distal sites to regulate cell-type spe- 
cific cytoskeletal genes using an enhancer-based strategy 
[26]. Second, since many non-coding RNAs (ncRNAs) are 
found in intergenic regions, it is possible that SRF regu- 
lates the expression of ncRNAs. To test this notion, I 
examined if there are annotated ncRNA loci within 
20 kb of the 1827 SRF intergenic binding sites. My in silico 
analysis showed that a significant portion (10.89%) of the 
intergenic peaks were associated with annotated ncRNAs 
(Figure 4D). A total of 316 unique ncRNAs (Additional 
file 11: Table S8) were associated with intergenic SRF peaks, 
indicating the possibility of regulation of these ncRNAs by 
SRF. A majority of the ncRNAs (84.49%, 267) identified are 
classified as long non-coding RNAs (IncRNAs), while 44 
are miRNAs (Figure 4E). In support of the possibility of 
SRF regulation of ncRNAs, the previously reported SRF tar- 
get miRNAs, miR-1 and miR-133, were on the miRNA list 
generated by my in silico analysis. 

UCSC Genome Browser images of SRF binding to 4 of 
its target genes {Srf, Vcl, c-fos and Egrl) are shown in 
Figure 5 A. Three representative MKLl peaks associated 
with Srf, Vcl and Cyr61 are shown in Figure 5B. Three 
representative ELK4 peaks on Srf, cfos and Egrl are shown 
in Figure 5C. 

MKL1 and ELK4 are dedicated to different transcription 
programs in macrophages 

To check the functional significance of DNA binding of 
SRF, MKLl and ELK4, I used the functional annotation 
tool, DAVID, to assess genes associated with binding peaks 
of SRF, MKLl and ELK4, respectively. The top 5 most 



relevant GO terms of biological process for each target 
gene list are presented in Figure 6A (full results in 
Additional file 12: Table S9, Parts I, II and III). SRF associ- 
ated genes were overrepresented with GO terms of "actin 
cytoskeleton organization" and "gene expression". The GO 
term "actin cytoskeleton organization" was also selectively 
enriched in MKLl associated genes. In contrast, "gene ex- 
pression" appeared as one of the most highly enriched 
terms for ELK4 associated genes. 

By using KEGG pathway analysis for the SRF, MKLl, 
and ELK4 bound genes, I found that MKLl and ELK4 
associated genes contribute to different pathways related 
to SRF. For example, SRF and MKLl associated genes 
showed overlapping and enriched pathways of "focal adhe- 
sion" and "adherens junction". In contrast, SRF associated 
genes were enriched in the "MAPK signaling pathway" 
with ELK4 associated genes. The KEGG term of "path- 
ways in cancer" was commonly seen in the analysis of all 
of the 3 groups of target genes (Figure 6B, full results in 
Additional file 13: Table SIO, Parts I, II and III). These 
functional differences can be illustrated by different target 
genes regulated by the two different cofactors. Actb is an 
example of a "focal adhesion" KEGG pathway molecule 
and Egrl is a representative "MAPK signaling pathway" 
molecule. Each of these genes are independently regulated 
by MKLl and ELK4 as indicated by genome browser im- 
ages. As shown in Figure 6C, Actb is an MKLl-dependent 
and ELK4-independent SRF target gene. In contrast, Egrl 
is bound by ELK4 and SRF, but not MKLl. 

I further confirmed the occupancy of MKLl and ELK4 
on different target genes by ChlP-qPCR and also com- 
pared the dynamic binding of MKLl and ELK4 before 
and after zymosan treatment. My experiments confirmed 
the abundant enrichment of MKLl on the promoters of 
5'// and Vcl, The observed enrichments were high under 
basal conditions and were further enhanced after zymo- 
san treatment (Figure 7A). This was in accordance with 
the increased expression of MKLl in response to zymo- 
san treatment in the BLRP-MKLl-BirA cell line. Con- 
ventional Chip analysis using primers flanking the peaks 
also confirmed ELK4 binding to the promoters of cfos 
and Egrl (Figure 7B). Although the mRNA expressions 
of these two target genes can be highly induced by zy- 
mosan treatment, ELK4 occupancy levels remained con- 
stant on these two promoters. 

Discussion 

In this study, I analyzed the role of SRF, and representa- 
tive members of its cofactor families in macrophages, in 
the context of zymosan stimulation. Using a knockdown 
strategy coupled with gene expression array analysis, I 
identified target genes whose expressions were dependent 
on the intact expression of SRF, MKLl/2 and ELK4, 
respectively. Chromatin immunoprecipitation combined 
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with high-throughput sequencing allowed the identifica- 
tion of direct target genes of SRF, MKLl and ELK4, Taken 
together, my results indicate that MKLl and ELK4 contri- 
bute to distinct aspects of SRF functions in macrophages. 

Here, I show that SRF mediates a response to zymosan 
stimulation in primary macrophages. The transcription 
of 5'// was increased by zymosan treatment in an early 
phase, similar to other highly characterized immediate 
early genes in the literature such as c-fos and Jun, Identi- 
fication of the SRF binding element, the CArG box, as 
one of the most highly enriched motifs in the promoters 
of zymosan induced genes suggests that SRF may play 
important roles in zymosan responses in macrophages. 



Using knockdown analysis combined with gene profiling, 
I identified different target gene sets for SRF and mem- 
bers of two well-known families of SRF accessory cofac- 
tors: the myocardin family members, MKLl/2, and the 
TCF family member, ELK4. With regard to the effects of 
knocking down SRF, I found that different gene sets with 
macrophage specific functions were targets of SRF, includ- 
ing genes related to migration, phagocytosis, adhesion and 
responses to stimuli. In agreement with previous studies 
performed in other cell types, MKLl/2 regulates the 
expression of many genes related to actin cytoskeleton 
organization in macrophages. These genes showed high 
basal expression and little to no induction in response to 
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Figure 7 Validation of MKL1 and ELK4 recruitment to target gene promoters following zymosan treatment. (A) MKLl ChIP assay was 
performed in BLRP-MKLl-BirA stable cells under basal conditions and after zymosan treatment for 1 h. Immunoprecipitated DNA was analyzed by 
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stable cells under basal conditions and after zymosan treatment for 1 h. Rabbit IgG was used as a negative control. Immunoprecipitated DNA was 
analyzed by qPCR using primers flanking the proximal promoter regions of c-fos and Egrl. 



zymosan treatment. In contrast, I found that ELK4 target 
genes are overrepresented in the GO term of "response to 
stress". I validated some representative ELK4 target genes 
and found that the expressions of these genes were highly 
induced by zymosan treatment and that their inductions 
were heavily dependent on intact ELK4 expression. Taken 
together, its proposed that MKLl/2 and ELK4 contribute 
to different functions in the quiescent state compared to 
the context of an inflammatory response. Moreover, the 
microarray data revealed that both MKL 1/2 and ELK4 have 
some SRF independent functions and that these functions 
may involve other novel cofactors. 

One of the limitations of microarray analysis is that it 
is difficult to differentiate the direct target genes from 
the secondary ones. For example, Srf expression itself 
can be affected by knocking down MKL 1/2 or ELK4. I 
also noticed that some ELK4 dependent target genes 
showed decreased expression with MKLl/2 loss, which 
is probably a secondary effect of SRF expression down- 
regulation. Therefore, I performed ChlP-seq analysis to 
further identify direct target genes of SRF and its cofac- 
tors and to explore novel partners of the three transcrip- 
tion factors. 

Motif analysis of the SRF bound peaks identified in my 
ChlP-seq experiment revealed a set of transcription fac- 
tors that could be cofactors of SRF besides ELK4. Some 



of these factors have already been reported as cofactors of 
SRF. For example, SRF acts in conjunction with TEADl 
for cardiac muscle gene regulation [38,42], the interaction 
of SRF and NFYA has been verified in vitro [39], and the 
interaction of SRF and NFkB has been observed on the 
IL-2R enhancer in T lymphocytes [40]. In addition, the 
NFkB subunit, p65 (also known as RELA), synergizes with 
SRF to activate an SRE containing reporter [41]. My study 
also revealed several additional potential cofactors of SRF, 
including GABPA, Staf, CREBl, STATl and RREBl. An 
independent ChlP-seq study in Jurkat T cells found that a 
large fraction of SRF peaks occur within 100 bp of GABP 
peaks, suggesting that SRF may interact physically with 
GABP [45]. Although CREB has not been reported as 
directly interacting with SRF, CREB binding protein 
(GBP) has been found constitutively bound to the SRE 
on the c-fos promoter [46]. RREBl and SRF were also 
both predicted to interact with the same factor, MEF2 
[43] . These references provide some evidence for my hy- 
pothesis that SRF may interact with cofactors outside of 
the myocardin and TCF families. 

Our previous ChlP-seq study of SRF in primary mac- 
rophages indicated a critical role of SRF in promoting 
hematopoietic-cell-restricted gene expression through bind- 
ing at distal enhancers, while modulation of generally ex- 
pressed genes occurred at promoters [26]. A large fraction 
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of the SRF peaks in primary quiescent macrophages also 
showed up in my ChlP-seq experiment in the BLRP- 
MKLl-BirA cells. Moreover, I also found that SRF-binding 
sites are not restricted to promoter regions, but are prefer- 
entially located in intergenic regions of the mouse genome. 
Through mapping of the intergenic peaks to the mouse 
genome, I found that 10.89% of the peaks were associated 
with annotated ncRNAs, of which, IncRNAs constitute the 
major fraction. A list of miRNAs were also associated with 
intergenic SRF peaks, including the known SRF targets, 
miR-1 and miR-133 [33]. While MKLl has a genomic dis- 
tribution pattern similar to SRF, as previously reported [19], 
ELK4 showed a relatively larger fraction of occupancy on 
promoters. In particular, GABPA and NFYA motifs were 
found to be enriched on ELK4 bound gene promoters, 
suggesting a possible partnership of the two transcription 
factors with ELK4. 

In addition to validating the MKLl and ELK4 binding 
sites from my ChlP-seq experiment, I also assessed the 
binding of these factors in the context of zymosan stimula- 
tion. With regard to mRNA expression, MKLl dependent 
genes showed high basal expression and little to no induc- 
tion in response to zymosan treatment. By ChlP-qPCR, I 
found that MKLl binding to the promoters of its target 
genes was increased after zymosan treatment. Considering 
the higher expression of MKLl after zymosan treatment in 
the BLRP-MKLl-BirA stable cells, this enhanced occu- 
pancy may reflect an increased expression of BLRP-MKLl 
protein. In contrast, ELK4 target gene mRNA levels can be 
highly induced by zymosan treatment. However, in the 
basal state, ELK4 is already bound to target genes poised 
for zymosan induction. I speculate that additional events, 
such as post-translational modification of the SRF-cofactor 
complex, occur upon zymosan treatment to promote the 
fiinction of ELK4-SRF at these loci. Kasza, et al reported 
similar results in which ELKl and SRF were both found on 
the ZC3H12A promoter and that their binding was also 
not increased upon IL-1|3 stimulation [47]. These results 
indicate that transcription factor recruitment may not be a 
key regulatory event in response to inflammatory stimuli. 

Conclusions 

In summary, this is the first study to examine the func- 
tions of SRF and its key cofactors simultaneously on a 
genome-wide level through both ChlP-sequencing and 
transcriptome analysis. The overall view of the genome- 
wide locations of SRF, MKLl and ELK4 has been con- 
structed and the results suggest that SRF utflizes MKLl 
to fulfiU steady state ceUular functions, including cyto- 
skeletal organization, and utflizes ELK4 to facilitate acute 
responses to external infection. Previous studies of SRF 
largely focused on skeletal, cardiac and smooth muscle 
cefls and neurons. This project sheds light on the fiinc- 
tions of SRF, MKLl and ELK4 in macrophages, indicating 



their potential roles for both cellular homeostasis and 
stress responses in immune cefls. 

Methods 

Cell culture and transient transfection 

Thioglycoflate (BD biosciences) -eUcited, peritoneal mac- 
rophages were prepared as previously described [48] from 
6-8 week old, male, C57BL/6 mice (Harlan). For RNAi 
experiments, 0.75 miflion primary macrophages were 
transfected with nonspecific (NS) control or SMART- 
pool siRNAs (100 nM, Dharmacon) directed against Srf, 
Midi, Mkl2 and Elk4 using the Deliver X transfection 
reagent (Panomics) according to the manufacturer s in- 
structions. RAW 264.7 were transiently transfected with 
a 3xCArG luciferase reporter using Superfect reagent 
(Qiagen) and luciferase activities were measured using 
the luciferase assay system (Promega). Cells were treated 
with or without zymosan A (1 mg/mL, Sigma- Aldrich) as 
stated in each experiment. 

Generation of BLRP-MKL1-BirA cell lines 

To achieve rapid and high affinity purification of MKLl 
for ChlP-sequencing assay, I generated a BLRP (biotin 
ligase recognition peptide)-TEV-Flag-MKLl expression 
vector containing a puromycin resistance gene. BLRP is 
a high affinity substrate for the E. coli biotin ligase, BirA, 
which biotinylates a lysine residue within the peptide 
[49,50]. A RAW264.7 cefl line that stably expresses BirA 
under neomycin selection has been previously described 
[49]. The BLRP-MKLl expression vector was trans- 
fected into the BirA expressing RAW264.7 stable cefl 
line and stable clones expressing both BirA and BLRP- 
tagged MKLl were identified by puromycin and neomy- 
cin double selection. 

RNA isolation and Q-PCR 

Total RNA (isolated by RNeasy kit, Qiagen) was pre- 
pared from untreated or zymosan treated primary mac- 
rophages and digested with DNase I according to the 
manufacturer's instructions. cDNA was synthesized from 
1 (ig of total RNA using Superscript III (Invitrogen) reverse 
transcriptase according to the manufacturer's instructions. 
Quantitative real-time PGR (Q-PCR) analysis was per- 
formed using 1 (il of cDNA template, 50 ng of each gene- 
specific primer, and SYBR greenER master mix (Invitrogen) 
on a 7300 real-time PGR system (Applied Biosystems). Pri- 
mer sequences are avaflable upon request. 

Gene expression profiling 

For the Mkll/2 double knockdown samples, 0.5 (ig of puri- 
fied RNA per sample was labeled using the LRILAK PLUS, 
two color low RNA input Linear Amplification kit and 
hybridized to Whole Mouse Genome Microarray 4 x 44 K 
60 mer slides (Agflent) according to the manufacturers 
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instructions. Slides were scanned using the Agilent GZ505B 
Scanner. For Elk4 knockdown samples, 100 ng of purified 
RNA per sample was reverse transcribed using M-MLV 
(Takara Chemicals) and amplified by an in vitro transcrip- 
tion based method. The Cyanine 3 (Cy3) or Cyanine 5 
(Cy5) labeled cDNA was hybridized to the CapitalBio 
Mouse Genome Oligo Array (CapitalBio) for 16 h at 42°C. 
Arrays were scanned with a LuxScan™ 10 K-A confocal 
scanner (CapitalBio), and the data were extracted with 
GenePix pro4.0 software (Axon Instruments). 

Gene ontology (GO) analysis and motif analysis 

GO analysis was performed using the web-based func 
tional-annotation tool, DAVID (http://david.abcc.ncifcrf. 
gov/home.jsp). Gene ontology terms were considered 
significant if they had a p-value of less than 0.05. The 
oPOSSUM human/mouse single site analysis tool (http:// 
www.cisreg.ca/oPOSSUM/) was used to identify motifs of 
12-bps in length that were overrepresented in different 
groups of target sequences [51]. Motifs with a Z-score of 
more than 5 and a Fisher Score of less than 0.05 were 
regarded as significantly enriched. Terms presented here 
are sorted by their respective Fisher scores. Analysis of 
the overrepresented miRNAs upstream of target genes 
was performed using Ingenuity Pathway Analysis (Ingenuity 
Systems, http://www.ingenuity.com). 

Chromatin immunoprecipitation (ChIP) assay 

Chip assays against MKLl were performed using the 
BLRP-MKLl-BirA cell line with methods similar to those 
previously described in [50]. Briefly, cells were crosslinked 
with formaldehyde then chromatin was fragmented into 
200-300 bp pieces by sonication and subjected to high af- 
finity purification with Dynabeads® MyOne™ Streptavidin 
Tl (Life technologies). The beads were subjected to 
AcTEV Protease (Life Technologies) digestion to release 
the tagged MKLl proteins complexed with interacting 
genomic DNA, leaving non-specific proteins bound to 
the beads. The MKLl-DNA elution without TEV diges- 
tion was used as a background control. Elutions from 
TEV digestion were reverse-crosslinked by heating, and 
DNA was recovered and purified. 

Chip assays against SRF and ELK4 were performed in 
the BLRP-MKLl-BirA cell line with methods similar to 
those previously described in [26]. Immunoprecipitations 
were performed by incubating with 2.5 (ig of anti-SRF 
antibody (sc-335, Santa Cruz Biotechnology) or anti-ELK4 
antibody (sc- 13030, Santa Cruz Biotechnology) overnight 
at 4°C. Protein-DNA complexes were recovered using 
50 \A of blocked ImmunoPure Protein A Agarose (Pierce) 
and then washed, eluted, digested and purified. Rabbit 
preimmune serum was used as a control for non-specific 
binding. 



High-throughput ChIP sequencing (ChlP-seq) and 
data analysis 

Purified ChIP DNA was adapter ligated and PCR amplified 
according to the manufacturers instructions (Illumina). 
High-throughput sequencing was performed on an Illu- 
mina Genome Analyzer and short reads (first 23 to 25 bp) 
were mapped to the mouse reference genome. Tag counts 
for each experiment were normalized to 10^ specifically 
mapped tags. Sequences with tag counts of more than 20 
were regarded as binding peaks. To determine SRF peaks, 
SRF ChlP-seq data was compared to a background data set 
where SRF ChlP-seq was performed in SRF knockout pri- 
mary macrophages. Only those peaks with tag counts at 
least 4-fold more than the tag counts in SRF knockout cells 
were considered specific. Peaks were visualized by pre- 
paring custom tracks using the UCSC Genome Brower 
(http://www.genome.ucsc.edu/) based on the protocol 
previously described [26]. Peaks were assigned to specific 
genes by proximity to the nearest transcription start site. 
To search for associated ncRNAs, sequences +/-20 kb 
away from target peaks were extracted and mapped to an- 
notated ncRNAs based on UCSC mm9 and NCBI build37. 

Additional files 



Additional file 1: Figure 51. SRF functions in the early phase of 
zymosan induction in macrophages. 

Additional file 2: Table SI. Genes that are downregulated (Part I) 
or upregulated (Part II) more than 1.5 fold in SRF knockdown primary 
macrophages after zymosan treatment for 1 h. 

Additional file 3: Table S2. Part I: Gene Ontology results (Biological 
process annotations, p < 0.05) for SRF downregulated genes 
(downregulated more than 1.5 fold with SRF knockdown in zymosan 
treated macrophages). Part II: Gene Ontology results (Biological process 
annotations, p<0.05) for SRF upregulated genes (upregulated more than 
1.5 fold with SRF knockdown in zymosan treated macrophages). 

Additional file 4: Figure S2. Knockdown of MKLl/2 and ELK4 has an 
impact on zymosan treated macrophages. 

Additional file 5: Table S3. Genes that are downregulated more than 
1.5 fold in MKLl and MKL2 double knockdown primary macrophages 
after zymosan treatment for 1 h. 

Additional file 6: Table S4. Genes that are downregulated more than 
1.5 fold in Elk4 knockdown primary macrophages after zymosan 
treatment for 1 h. 

Additional file 7: Figure S3. Construction of BLRP-MKLl-BirA stable 
cell lines. 

Additional file 8: Table S5. SRF binding peaks in BLRP-MKLl-BirA stable 
cells (tag count more than 20). 

Additional file 9: Table S6. MKLl binding peaks in BLRP-MKLl-BirA 
stable cells (tag count more than 20). 

Additional file 10: Table S7. Elk4 binding peaks in BLRP-MKLl-BirA 
stable cells (tag count more than 20). 

Additional file 11: Table S8. The 316 ncRNAs associated with SRF 
intergenic peaks. 

Additional file 12: Table S9. Part I: Gene Ontology results (Biological 
process annotations, p < 0.05) for SRF peak associated genes (tag count > 20, 
in zymosan treated BLRP-MKLl-BirA cells). Part II: Gene Ontology results 
(Biological process annotations, p < 0.05) for MKLl peak associated genes 
(tag count > 20, in zymosan treated BLRP-MKLl-BirA cells). Part III: Gene 
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Ontology results (Biological process ar^r^otatior^s, p < 0.05) for ELK4 peak 
associated genes (tag count > 20, in zymosan treated BLRP-MKLl-BirA cells). 

Additional file 13: Table S10. Part I: KEGG pathway analysis full results 
(p < 0.05) for SRF peak associated genes (tag count > 20, in zymosan 
treated BLRP-MKLl-BirA cells). Part II: KEGG pathway analysis full results 
(p < 0.05) for MKLl peak associated genes (tag count > 20, in zymosan 
treated BLRP-MKLl-BirA cells). Part III: KEGG pathway analysis full results 
(p < 0.05) for ELK4 peak associated genes (tag count > 20, in zymosan 
treated BLRP-MKLl-BirA cells). 
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